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We investigate the peculiar velocity field due to long cos- 
mic strings in several cosmological models and analyse the 
influence of a nonscaling behaviour of the string network, 
which is expected in open cosmological models or models with 
a cosmological constant. It is shown that the deviation of 
the propability distribution of the peculiar velocity field from 
the normal distribution is only weak in all models. It is fur- 
ther argued that one can not necessarily obtain the parame- 
ter f3 = fij]' 6 jb from density and velocity fields, where fio is 
the density parameter and b the linear biasing parameter, if 
cosmic strings are responsible for structure formation in the 
universe. An explanation for this finding is given. 



I. INTRODUCTION 

To understand the origin and the formation of struc- 
ture in the universe is one of the most challenging prob- 
lems in modern cosmology. There are two competitive 
theories which try to explain the origin of the seeds. The 
first one is inflation, in which the universe undergoes 
an epoch of fast (inflationary) expansion, triggered by 
a scalar field, called the inflaton. Quantum fluctuations 
in the inflaton field are stretched on supcrhorizon scales 
which arc turned into matter fluctuations at later times. 
These matter fluctuations represent the seeds of the ob- 
served structure today In the other theory topolog- 
ical defects are responsible for structure formation. For 
example, cosmic strings might be produced in a phase 
transition in the very early universe. If they are heavy 
enough, they influence the cosmological fluid and could 
seed the structure in the universe (see |2|J^] for reviews 
and references). 

There are several ways to test these theories. For 
example, they make different predictions for the angu- 
lar fluctuation spectrum of the CMBR on small scales. 
The predicted spectra can then be compared with the 
data, for example of future projects such as MAP and 
PLANCK. Related to this is the comparison of the pre- 
dicted matter power spectrum with the data. Cosmic 
strings could further be tested with "astrophysical" tests, 
such as the expected gravitational radiation background 
from strings etc. It's interesting to note, that all these 
tests gave consistent results for the string parameter /z, 
the mass per length on the string 



Another possibility for testing structure formation the- 
ories was proposed by Q]. The propability distribution 
of the peculiar velocity field should be different in infla- 
tionary models and models with topological defects such 
as cosmic strings. However, as emphasized by |^||, the 
propability distribution of the peculiar velocity field in 
cosmic string theories is Gaussian to high accuracy. This 
conclusion was based on the assumption that the string 
network reaches a scaling behaviour. 

It was shown by several authors that the scaling so- 
lution is only expected in the Einstein-de Sitter model 
[Mpl; m °P en models, in flat models with a cosmological 
constant and in closed (loitering) models the behaviour 
of the network is different from scaling. Further, the 
transition to the matter scaling behaviour is much longer 
than previously estimated These possible sources of 
deviation from the string scaling solution should have in- 
teresting consequences. A first step was done in |To| ] and 
JO} . Whereas in Jh} it was shown, that only a drastic 
departure from scaling could solve the problems of struc- 
ture formation which cosmic strings, |TT| have shown that 
there is only weak dependence of the density parameter 
f^o m open and flat models with cosmological constant on 
the normalisation of G/i from COBE data. It might be, 
that the cosmic string scenario is successfull also in the 
absence of the scaling behaviour, i.e. that it works well in 
open models or in models with a cosmological constant. 

In this paper we investigate the influence of cosmic 
strings on the peculiar velocity field. Earlier investiga- 
tions of the peculiar velocity field concentrated on the 
spectrum of the field, i.e. its dependence on the length 
scale L, see e.g. P, ^2| , ^3[ . We are interested in the ef- 
fects of a departure from the scaling behaviour. We use 
an approximation, first introduced by Q to calculate the 
effects of many strings. 

The paper is organized as follows: In section II we 
discuss the influence of cosmic strings on the peculiar ve- 
locity field. Our calculations of the string network are 
based on the calculations by [pj. Our results for the pe- 
culiar velocities are presented in section III. In section 

IV we argue that if cosmic strings seed the structure in 
the universe then the peculiar velocity field and the den- 
sity field is correlated but one cannot obtain information 
on the parameter (3 = fl^/h, where fio is the density 
parameter and b is a linear bias parameter. In section 

V we summarize our results and give some conclusions. 
Throughout the paper we set c = I . 



1 



II. THE PECULIAR VELOCITY FIELD DUE TO 
LONG COSMIC STRINGS 

The space-time of a straight cosmic string is similar to 
the Minkowski space-time, except for a deficit angle A</>, 
given by 



A(f> = 8nGp:j s v s 



(2.1) 



As a result, the matter gets a kick towards the plane 
swept out by the string (v s is the string velocity and "f s 
is the Lorentz factor). The velocity kick due to a wiggly 
string is given by 



u s = AirGfi^ s v s f = 3.8(Gp) 6 (j s v s )f km/s, (2.2) 



with 



1 
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(2.3) 



The term / corresponds to the small scale structure on 
the string, where \x is the effective mass per unit length 
on the string and T is the string tension. 



A. The Zeldovich approximation 

To calculate the peculiar velocity field, we use the Zel- 
dovich approximation, in which the physical coordinates 
of a particle are written by 



r(x,t) = a(t) [x + V(x,t)] , 



(2.4) 



where a(t) is the scale factor, x is the comoving coordi- 
nate and ip is the displacement vector due to inhomogen- 
ities in the cosmic fluid, i.e. cosmic strings in our context 
mM . The equation of motion is given by Newtons law 



(2.5) 



The gravitational field is connected with the matter dis- 
tribution (Poisson equation), which can be obtained from 
linearising Einstcins field equation: 



V 2 .$ = 4 7 rG(p 6 + <5 j o) + Ac 2 



(2.6) 



In this equation pt, is the matter density, Sp the matter 
density fluctuation and A the cosmological constant. To 
first order one obtains 



6 =P_Pb = _ Vx .^ (x t ) 
Pb 



(2.7) 



Where p is the total matter density. This leads to 

V 2 $ = 47rGp h (l - V x -i/;)+Ac 2 (2.8) 
Integration of this equation and substitution of 



r = — r + 2aip + dip 
a 



and the second Friedmann equation 



o 4ttG Ac 2 
a 3 3 



(2.10) 



leads to the evolution equation for the displacement ip 



t/j + 2- ip-AnGpb ip = 0. 



(2.11) 



For our purposes we have to calculate the pecu liar veloc- 
ity field, which can be obtained from eq. ( [2.4|) : 



v pec 



aip. 



(2.12) 



The effect of the cosmological constant on the evolution 
of a density perturbation is only due to the effect of A on 
the evolution of the scale factor a. 



B. The influence of cosmic strings 

We use an analytical approximation (the so called mul- 
tiple impulse approximation) , first introduced by Vachas- 
pati |J, which also was successfully applied to calculate 
the CMBR anisotropies pH] . We devide the time inter- 
val from t eq (at which structures starts to form) to to 
in N steps with ti + \ = 2ti. Between ti and t; + i the 
strings intercommute, form loops etc. so that (approxi- 
mately) at tj_|_i the "new ordered state" of the network 
is uncorrelated with the "old state" at ti. Again, at this 
time the network influences the matter within the horizon 
(due to scalar field radiation scales larger than the hori- 
zon are not affected). This is not true in vacuum domi- 
nated epochs. In this case the velocity of the strings de- 
creases and therefore the propability of string interaction 
decreases. This means that the new state is correlated 
with the old one. However, in such epochs the number 
of strings within the Hubble horizon decreases rapidely 
and therefore our results aren't changed significantly (see 
below). 

At ti = t eq each string within the horizon gives the 
matter a kick in the direction of the surface swapped out 
by the string: 



^s^-i' 



(2.13) 



where k^i is a (random) unit vector in direction of the 
string i. The resulting peculiar velocity from all strings 
at t\ is 



vi 



= ^ u ski,i, 



(2.14) 



The sum is now taken over the number n 9y i of all strings 
within the horizon at t\. This peculiar velocity field 
grows between t\ a nd t^ by a factor Aitx^t^) via eq. 
( 2 - 9 ) ( ETill ) and eq. (§T|): 
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A(ti,t f ) 



a(tf)m f )\ 



(2.15) 



a(ti)\ip(ti)\ 

At the time t 2 the peculiar velocity field is given by 

"3,2 

2 n 9 ,i 

2£^(*a»*i)kij ■ ( 2 - 16 ) 
_i=i 3=1 

Here we have used that A(ti,tj)A{tj,t{) — A(ti,t{) and 
that the velocity of the strings is the same in every epoch. 
For our purposes this is a good approximation, because 
when the strings slow down the number of strings within 
the horizon also decreases. 

Iteration leads to the peculiar velocity field at the 
present time on a scale L: 



ML) 



Y^^2A(t ,ti)k4 

i=l 3=1 



(2.17) 



In this equation Nl is the number of Hubble time steps 
during which a volume of comoving size L 3 experiences 
string impulses, n g j is the number of strings within the 
horizon at the time tj. We assume that the vectors kjj 
are random, that is 



< kjj • k; /m >= SuS 



(2.18) 



III. RESULTS 

We calculate the peculiar velocity field for four repre- 
sentative models, shown in Tab. 1. For the Einstein-de 
Sitter (E-dS) model we discuss the influence of the long 
transition between the radiation and matter scaling solu- 
tion [fl. As a test, we include the case for an ideal scaling 
in the E-dS model. The propability distributions of the 
peculiar velocity fields at a scale L(t eq ) for the models 
are shown in Fig. 1-5. Each plot was obtained after 
50000 realisations. For an exact scaling behaviour in the 
E-dS model this distribution was shown to be Gaussian 
||||. We obtain the same result (see Fig.l). In the case 
of the long transition between the radiation and matter 
scaling behaviour the distribution remains nearly Gaus- 
sian. However, the distribution becomes broader and the 
peculiar velocity increases (Fig. 2). The Gaussian char- 
acter of the propability distribution can be found in the 
other models, too. There is only a slight deviation at 
large and small velocities. 

In the models, we obtain a peculiar velocity at a scale 
corresponding to L(t eq ) given by (the length scales are 
calculated with Hq = 100km/(s-Mpc)): 

v pec {L eq ps 70 Mpc) closc d (3.1) 
= (460 ± 200){Gfi) 6 { ls v s )f km/s, 

(L eq » 1 Mpc) Eds , ni (3.2) 
= (1740 ± 760) (G/i) 6 ( 7s u s )/ km/s, 



From these equations we calculate the RMS velocity nu- 
merically on a scale L(t eq ). On scales smaller than L{t eq ) 
the peculiar velocity field depends only weakly on L 
whereas on scales larger than L(t eq ) the predicted ve- 
locity field scales as L^ 1 [[l2[. 



C. Network parameter 

We use the calculation from for the statistical prop- 
erties of the string network. We set the number of strings 
within the horizon H~ x by 



n a = l + {£>H)- 



(2.19) 



where £ is the characteristic length scale of the string 
network, defined by 



(2.20) 



where is the density of the long strings. In the radia- 
tion dominated epoch n s is about 10, in the (late) matter 
dominated epoch (with scaling) this number is about 3. 
In more general cosmological models this number is a 
function of ti me |8|| . Later we will discuss the influence 
of the ansatz (|2.19|). 



(L eq ps 1 Mpc) EdSi id 
- (1240±570)(G/i)6(7^ s )/ km/s, 

V pec {L eq P3 10 Mpc) A ,flat 

= (280 ± 120)(G/*) 6 (7.«.)/ km/s, 



(L^ « 10 Mpc) 

open 

(80±35)(G/i) 8 (7.i>.)/ Ws- 



(3.3) 



(3.4) 



(3.5) 



The length scale corresponding to the time t eq is set to 
be O.Lff(^)- 1 [§: 



: e9 pa 1.1-VMpc 

"0 



(3.6) 



Note, that in || the length scale was set to be 0.7t e q- 
Therefore, we assume a somewhat pessimistic view when 
strings could significantly influence the volume at t eq . 
In our picture, the volume must be within the typical 
length scale between all strings. However, the volume is 
influenced by the strings outside the volume and there- 
fore we somewhat underestimate the peculiar velocities. 
However, this can be taken into account with including a 
para met er £, which modifies our ansatz ( 2.19| ) (see below, 
eq. (ph). 



3 



Within this length the peculiar velocity remains nearly 
constant, because a smaller length corresponds to times 
t < t eq in which perturbation grow only weakly. This 
would imply that for the closed model we would expect 
nearly constant bulk flows on scales smaller than 70 Mpc, 
which is indeed observed. The situation in the other 
models is not so clear, because for scales larger that L eq 
the velocity decreases as L increases (v oc L^ 1 ). 

It is interesting to note that in all models the standard 
deviation is related to the mean value by 



a « 0.45u n 



(3.7) 



For our calculations we have used the ansatz ( 2.19| ) 
for the number of strings within the horizon. Although 
this should be a good approximation we could set n s = 
C(l + (£ ■ H) )■ The frequency distribution remains 
Gaussian, however, now the RMS velocity and the stan- 
dard deviation is given by 



VpecX = \/C v pecX=l 



and 



o- C =i- 



(3.8) 



(3.9) 



Here, v pec x=i and is given by eq. (3.1-3.5) for the 
cosmological models. The peculiar velocities therefore 
depend on the parameter 



= VCV(>{Vsls)f- 



(3.10) 



To conclude, the fluctuation of number of strings 
doesn't change the shape of the propability distribution 
of the peculiar velocities and the amplitude depends on 
the same set of parameters (3.11) as in the case for a 
ideal scaling behaviour of the string network. However, 
the effective number of strings and the maximum length 
on which coherent bulk flows are expected, depends on 
the cosmological parameters £l m ,o, an d Hq. 



IV. MATTER DISTRIBUTION, BULK FLOWS 
AND BIASING 

The results presented in the last section imply that 
the parameter j3 — fl®f /b could not be obtained from 
velocity-density reconstruction methods such as PO- 
TENT. To see this, we remember that the fundamental 
equation, on which these kinds of reconstruction methods 
are based, is given by [ p^| , p7| 



V • v = -(3H5. 



(4.1) 



Here, 8 is the density fluctuation. On the other side, the 
continuity equation holds: 



= -8. 



(4.2) 



In fact, in line ar approximation eq. (4.1) can be obtained 
from eq. (4.2). The important point is that if the ratio 
8/8 is independent of space, an arbitrary application of 
eq. (4.1) can lead to an under- or overestimation of (3 
if one applies eq. ([4.l[ ) arbitrarily to the data. This was 
shown by jl^] in the context of the explosion scenario. To 
demonstrate this point we repeat their short analytical 
example: 

Let us consider an empty universe with Qq = 0, filled 
with massless particles. At some time ti, the matter gets 
a kick due to a cosmic string (in the paper by Babul et al., 
the case of the explosion scenario was considered, but in 
the case of cosmic strings the analysis is identical). The 
linear Euler equation reads: 



~dt 



2H(t)v = v stri ng<5(£ - U). 



The density contrast evolves according to 
d 2 8 , . 88 



(4.3) 



(4.4) 



This equation can be solved with the boundary conti- 
tions at the time ti, which are <5(x, ti) = and 8(x, ti) = 
C(x,t } ) = — V x • v str ing(x, ti). The solution of the equa- 
tion (L4) is given by (a(t;) = 1) 



<5(x,t) = C(x,ti)*i(a(*) - l)Mt), 

t<y(x,*) = c(x, 



(4.5) 
(4.6) 



One can see, that the ratio 8/8 is independent of space. 
We can use the continuity equation (^^) and equation 
(0) to get 



PeS = 



1 



a[t) - 1' 



(4.7) 



Alth ough the true value is zero, an observer, applying 
(11) to the data, will get a value that is different from 
the true one. Only at late times /3 c ff will approach the 
value [3 — 0. 

We expect similar results for the models in table 1. 
The cosmological model changes the time dependence of 
8 and 8, but in general, if there were only one velocity 
kick on the matter, the ratio 8/8 would be independent 
of space. 

We have seen, that in cosmic string models of struc- 
ture formation the observed peculiar velocity field is a 
vector sum of many contributions of strings. We write 
the general solutions of 8 and 8, which arise from one 
kick as 



SijfaU) = C(^U)i,jB(t,ti), 
Sijfati) = ((x,ti)i^C{t,ti). 



(4.8) 
(4.9) 



The fields <5 and 8, which arise from all strings in the past 
are given by 
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i=i i=i 

AT n B ,i 

«(x,t) = ^;£^(x,t). 

i=l 3=1 

If we now replace Qj by 8i j/B, we can write 



(4.10) 
(4.11) 



N n g ,i 



s = X! ^2 kj ( x ' il 



i=l 3 = 1 



C(t,U) 
B(t,U)' 



(4.12) 



If we compare this with eq. (4.8) we see that there is a 
relation between the velocity field and the density field, 
but not of the same form as in eq. fl4.1| ). For cosmic 
strings we could write 



JV Kg,- 

z=l 3=1 



C(t,U) 
B(t,ti)' 



(4.13) 



If we compare eq. (4.13) with (4.1), we found that 

1 Hi, 3 ^i,fii,3 /BiJ 



fleff = 



(4.14) 



j "i,3 



where B and C are the solutions of the perturbation equa- 
tions. This equation can be interpreted as follows: each 
kick g ives an effective /3 ff,i,j = Cij/Bij. Therefore, eq. 
(4.14) is the weighted mean of the P e s,i,j- 



(3eS 



1 Ei, 7 - kiPeft, 



i,3 



h y, K , Si,j 



(4.15) 



This represents a measure of the departure of an exact 
relation between the velocity and density field. 



V. DISCUSSION & CONCLUSIONS 

In this paper we have considered the properties of the 
peculiar velocity field of galaxies in structure formation 
theories with cosmic strings. We considered the fact that 
the string network might not have developed a scaling 
behaviour (as is the case in open models or models with 
a cosmological constant) and showed that the propabil- 
ity distribution of the peculiar velocities is nearly Gaus- 
sian. The RMS peculiar velocity depends on the (effec- 
tive) number of strings within the horizon and on the 
string parameter G\x. The length, within the peculiar 
velocity is nearly independent of the scale, depends on 
the cosmological parameter f2o; and Hq. 

Open models have more problems with the amplitude 
of the peculiar velocity field. Only an unplausible high 
value of / could solve the problems (/ ~ 5). The situa- 
tion might be better in flat models with a cosmological 
contant. However, on scales larger than L eq the peculiar 
velocity drops with L , that is on a scale of about 60 



Mpc we expect in model 4 peculiar velocities of 50-100 
km/s. This is not in agreement with the observed value 
of 350-450 km/s. The situation is very good in the closed 
model. On scales smaller than L eo the velocity increases 
only weakly as L decreases |5| JT2] , [l3]] and remains nearly 
constant at 460 km/s (for Gfi = 10~ 6 and f"f s v s ~ 1) up 
to scales of about 200 Mpc. 

The results imply that if cosmic strings seed the struc- 
ture in the universe, one can not necessarily obtain the 
density parameter from the data. Comparison of density 
fields and velocity fields lead to an effective value, which 
is a measure of the deviation of an exact relation between 
the velocity and density field. 

Further work should be done on structure formation 
with cosmic strings in order to investigate the effects on 
a non-scaling behaviour of the cosmic string network. 
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Model 


K 


Ho 


A 


Ho 


N eq 


1 


+1 


0.014 


1.08 


90 


13 


2 





1.0 


0.0 


60 


20 


3 


-1 


0.1 


0.0 


60 


14 


4 





0.1 


0.9 


60 


15 



TABLE I. The four representative cosmological models. K 
is the curvature parameter, Oo is the matter density parame- 
ter, Ao is the cosmological term, Ho is the Hubble parameter 
(in km s _1 Mpc _1 ). N eq is the number of Hubble steps be- 
tween t eq and to- 
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v/us 



FIG. 1. The propability distribution of the peculiar veloc- 
ity field at a scale corresponding to t eq in the E-dS model 
with ideal scaling. The solid curve is the normal distribution. 
The mean value of v pec /u s is given by 325 and the standard 
deviation by 150. 
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v/us 

FIG. 2. The propability distribution of the peculiar veloc- 
ity field at a scale corresponding to t eq in the E-dS model with 
non-ideal scaling. The solid curve is the normal distribution. 
The mean value of v pec /u s is given by 460 and the standard 
deviation by 200. 
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v/us 

FIG. 3. The propability distribution of the peculiar veloc- 
ity field at a scale corresponding to t eq in the closed model. 
The solid curve is the normal distribution. The mean value 
of Vp ec /u s is given by 123 and the standard deviation by 55. 
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FIG. 4. The propability distribution of the peculiar veloc- 
ity field at a scale corresponding to t eq in the flat model with 
cosmological constant. The solid curve is the normal distri- 
bution. The mean value of v pec /u s is given by 73 and the 
standard deviation by 33. 
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FIG. 5. The propability distribution of the peculiar veloc- 
ity field at a scale corresponding to t eq in the open model. 
The solid curve is the normal distribution. The mean value 
of Vp ec /u s is given by 22 and the standard deviation by 9.8. 
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